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PREFACE 
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element  method  of  analysis  for  the  analysis  of  concrete  gravity  dams.  This 
Phase  II  report  will  address  only  the  dynamic  analysis  of  the  gravity  dam  and 
the  need  for  including  a  foundation.  The  Phase  la  report  addressed  the  static 
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tion  in  the  static  analysis  of  concrete  gravity  dams.  Other  reports  will  ad¬ 
dress  guidance  for  other  phases  of  finite  element  analysis.  The  work  was 
sponsored  under  funds  provided  to  the  US  Army  Engineer  Waterways  Experiment 
Station  (WES)  by  the  Civil  Works  Directorate,  Office,  Chief  of  Engineers 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

feet 

inches 

kips  per  foot 

kips  per  square  foot 

pounds 

pounds  per  cubic  foot 

pounds  per  foot 
pounds  per  square  foot 
pounds  per  square  inch 
slugs 


_ 1Z _ 

0,3048 

2.54 

1355. 818 
47.88026 
4.448222 
16.01846 

14.5939 

47.88026 

6.894757 

14.5939 


_ To  Obtain _ 

metres 

centimetres 

newton-metres 

kilopascals 

newtons 

kilograms  per 
cubic  metre 

newtons  per  metre 

pascals 

kilopascals 

kilograms 
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THE  RESPONSE-SPECTRUM  DYNAMIC  ANALYSIS  OF  GRAVITY  DAMS 
USTNG  THE  FINITE  ELEMENT  METHOD 


Phase  II 


PART  I :  INTRODUCTION 


Background 


1.  This  report  has  been  prepared  as  part  of  the  ongoing  effort  by  the 
Computer-Aided  Structural  Engineering  (CASE)  Committee  on  finite  element  anal¬ 
ysis.  It  is  part  of  a  Corps-wide  project  to  provide  guidance  for  the  use  of 
finite  element  analysis. 


Objective 

... 

2.  The  primary  objective  of  this  study  is  to  give  direction  in  perform¬ 
ing  a  response-spectrum  dynamic  finite  element  analysis  of  a  gravity  dam  using 
a  general-purpose  computer  program.  This  is  one  of  several  methods  currently 
being  used  within  the  Corps  of  Engineers  (CE)  for  performing  a  dynamic  analy¬ 
sis  on  gravity  dams  based  on  guidance  provided  in  ETL  1110-2-303  (Department 
of  the  Army  1985)  .  Other  procedures  and  computer  programs  commonly  used  are 
'  outlined  below  with  the  differences  between  them  discussed.  Table  1  contains 

*\a  summary  of  the  procedures  and  computer  programs. 

» 

a.  The  seismic  coefficient  method  as  presented  in  EM  1110-2-2200 
(Department  of  the  Army  1958)  is  in  reality  a  static  analysis 
with  static  forces  representing  inertial  and  hydrodynamic  loads. 
The  dam  is  analyzed  using  elementary  beam  theory  and  assumes  a 
fixed  foundation  and  incompressible  reservoir.  The  hydrodynamic 
pressure  distribution  increases  with  depth.  This  procedure  is 
good  for  horizontal  motions  only.  The  procedure  is  currently 
used  only  to  evaluate  overturning  and  sliding  stability. 

b.  The  "simplified"  method  for  earthquake  analysis  of  gravity  dams 
was  developed  by  A.  K.  Chopra  (1978).  The  analysis  of  the  dam 
is  similar  to  that  of  the  seismic  coefficient  method  in  that  the 
inertial  and  hydrodynamic  loads  are  represented  as  static  forces 
and  the  dam  Is  analyzed  using  elementary  beam  theory.  The  sim¬ 
plified  method  also  assumed  a  rigid  foundation.  The  difference 
between  the  seismic  coefficient  method  and  the  1978  version  of 
the  simplified  method  lies  in  the  way  that  the  inertial  and 
hydrodynamic  loads  are  determined.  In  the  simplified  method 
these  loads  are  based  on  the  fundamental  mode  of  vibration  of 
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Table  1 

Currently  Used  Procedures  Within  the  Corps  of  Engineers  for  Dynami 

Analysis  of  Gravity  Dams 
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the  dam,  with  and  without  reservoir.  In  this  method  the  reser¬ 
voir  is  assumed  to  be  compressible.  As  with  the  seismic  coeffi¬ 
cient  method,  only  horizontal  motions  are  considered.  In  1986, 
Fenves  and  Chopra  revised  the  simplified  method  to  include  the 
effects  of  an  elastic  foundation  (by  using  a  period  lengthening 
ratio  and  added  damping)  and  the  effects  of  sediments  in  the 
reservoir . 

£.  In  1986,  Cole  and  Cheek  (Technical  Report  SL-86-44,  Department 
of  the  Army)  combined  the  simplified  method  (1978  version)  with 
a  finite  element  model  of  the  dam  (in  lieu  of  using  elementary 
beam  theory) .  At  present  the  procedure  assumes  a  rigid  founda¬ 
tion  but  plans  are  to  include  the  effects  of  a  linear  elastic 
foundation  modeled  by  a  finite  element  grid. 

d.  EADHt  (a  computer  program  for  Earthquake  Analysis  of  Gravity 
Dams  Including  Hydrodynamic  Interaction)  is  a  finite  element 
program  that  uses  a  time-history  (but  in  the  frequency  domain) 
method  of  dynamic  analysis.  The  program  assumes  a  rigid  founda¬ 
tion  and  models  the  reservoir  as  a  comnressible  fluid.  The 
hydrodynamic  effect  is  represented  by  applying  added  mass  to  the 
face  of  the  dam.  The  added  mass  is  computed  using  methods  simi¬ 
lar  to  those  used  in  the  simplified  method  discussed  earlier. 

e.  The  EAGD-84  (Earthquake  Analysis  of  Gravity  Dams  Including 
Hydrodynamic  and  Foundation  Interaction)  program  is  verv  similar 
to  the  EADHI  program  except  that  the  foundation  is  modeled  as  a 
viscoelastic  halfplane. 

f.  Several  general-purpose  programs  are  available  that  will  perform 
either  a  2-D  or  3-D  finite  element  time-history  or  response- 
spectrum  dynamic  analysis.  Both  dam  and  foundation  can  be  mod¬ 
eled  by  a  linear  elastic  finite  element  grid  or  the  dam  can  be 
modeled  on  a  fixed  foundation.  The  influence  of  the  reservoir 
on  the  dam  is  typically  modeled  using  the  concept  of  added  mass. 
Westergaard ' s  added  mass  method  is  commonly  used  which  assumes  a 
rigid  dam  and  incompressible  fluid. 

3.  This  study  is  a  continuation  of  the  Phase  la  study  which  familiar¬ 
ized  the  beginning  finite  element  analyst  with  the  steps  necessary  to  perform 
a  static  fir.itc  element  analysis  cf  a  gravitv  dam  (Will  1987).  Many  of  the 
steps  in  performing  this  dynamic  analysis  could  also  apply  to  other  Corps 
structures.  The  beginning  engineer  should  be  able  to  develop  an  understanding 
of  the  steps  necessary  to  perform  a  dynamic  analysis  by  carefully  following 
this  study  and  using  supplemental  material. 

4.  Fallowing  is  a  list  of  the  necessary  steps  in  performing  a  dynamic 
finite  element  analysis.  The  steps  are  very  similar  to  those  for  a  static 
analysis  as  presented  in  the  Phase  la  report  (Will  1987). 

a.  Select  a  finite  element  computer  program  currently  in  use  by  the 
Corps  or  in  widespread  use  by  private  engineering  firms  and  sup¬ 
ported  by  a  vendor  with  the  desired  analysis  capabil ities. 
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b.  .elect  a  simple  problem  as  close  as  possible  in  overall  geome¬ 
try,  material  properties,  boundary  conditions,  and  loading  con¬ 
ditions  as  the  real  structure  to  be  analvzed.  This  structure 
should  have  closed-form,  experimental ,  cr  other  analytical  so¬ 
lution  results  available. 

£.  Select  the  type  of  analysis  (i.e.  response-spectrum,  time- 
historv)  to  be  performed  to  obtain  the  desired  results. 

d.  Select  the  finite  element  types  to  be  used  in  the  analysis  from 
the  library  of  elements  available  in  the  program  chosen  in 
Step  a. 

e.  Develop  and  analyze  finite  element  models  of  the  simplified 
structure  and  compare  results,  such  as  deflections  and  stresses, 
with  the  closed-form  results. 

_£.  Develop  modeling  guidelines  for  both  the  grid  and  loadings  from 
the  results  of  Step  e,  which  mav  be  extended  to  the  real 
structure . 

g.  Prepare  a  finite  element  model  of  the  real  structure  and  perform 
an  analysis. 

h.  Ask  the  following  question:  Ts  the  solution  accurate?  If  the 
answer  is  no,  refine  and  reanalyze  until  the  answer  is  yes. 

s.  Before  actually  performing  the  analysis,  further  detailed  discussion 
of  these  steps  is  warranted  to  understand  their  necessitv: 

a.  In  Step  a,  the  key  concept  is  that  the  finite  element  program 
should  be  currently  used  by  the  Corps  or  other  engineering  firms 
and  supported  by  a  vendor.  There  are  numerous  finite  element 
programs  available  today,  therefore  care  must  be  taken  in  the 
selection  process.  While  factors  such  as  ease  of  use,  func¬ 
tional  capabilities,  and  price  are  extremely  important,  an  over¬ 
riding  consideration  is  the  use  of  the  program  within  the  Corps 
or  other  engineering  firms  and  support  by  the  vendor.  An  ideal 
situation  is  to  find  a  program  that  is  easy  to  use,  has  the  nec¬ 
essary  functional  capabilities,  is  reasonably  priced,  and  is 
currently  being  used  by  someone  within  the  engineer's  group  and 
is  supported  by  the  vendor. 

b.  The  motivation  for  Steps  a-e  is  to  provide  an  opportunity  for 
the  engineer  to  build  confidence  in  the  use  of  the  program,  fi¬ 
nite  element  modeling  techniques,  and  to  develop  an  understand¬ 
ing  of  the  convergence  criteria.  Another  important  reason  for 
these  seeps  is  to  provide  the  engineer  with  an  understanding  of 
the  type,  quantity,  and  quality  of  finite  element  results.  A 
much  too  common  occurrence  is  for  the  engineer  to  devote  an 
enormous  amount  of  time  in  developing  the  finite  element  model 
and  after  results  have  been  obtained,  too  little  time  is  devoted 
to  the  interpretation  of  these  results,  i.e.,  the  accuracy  of 
the  results,  or  their  actual  usage  in  the  design  process. 

£.  From  the  analysis  performed  In  Steps  b-e ,  the  engineer  must  then 
extrapolate  the  information  gained  from  the  modeling  of  the  sim¬ 
ple  structure  to  the  modeling  of  the  real  structure.  Guidelines 


such  as  the  number  of  subdivisions  of  the  mesh  in  the  horizontal 
and  vertical  directions  may  be  developed  for  use  in  the  first 
model  of  the  real  structure.  Extreme  care  should  also  be  used 
when  deciding  how  to  handle  the  mass  distribution  of  the  struc¬ 
ture  and  the  necessity  to  include  any  additional  mass,  such  as 
foundation  and  adjacent  water. 

d.  In  Step  g,  the  real  structure  is  modeled,  analyzed,  and  the  re¬ 
sults  are  interpreted.  This  leads  to  the  crucial  questions  in 
the  analysis:  Ts  the  solution  accurate  within  an  error  criteria 
developed  by  the  engineer?  How  much  error  is  there?  These  are 
the  most  difficult  and  crucial  questions  in  the  entire  process. 
In  many  instances,  the  only  correct  way  to  answer  these  ques¬ 
tions  is  to  refine  the  model,  reanalyze,  and  compare  solutions. 
The  following  question  should  then  be  asked:  Have  the  results 
changed  significantly  due  to  the  refinement?  If  not,  an  approx¬ 
imate  solution  has  converged  and  the  engineer  must  determine  if 
the  results  make  physical  sense.  If  the  results  have  changed 
significantly,  other  models  may  be  required  and  comparisons  re¬ 
peated  until  convergence  is  satisfied.  The  engineer  must  remem¬ 
ber  that  the  finite  element  method  is  an  approximate  solution 
technique . 

e.  Also  the  question  of  whether  to  perform  a  time-history  or 
response-spectrum  analysis  should  be  considered.  In  either  case 
the  dynamic  loadings  and  mass  distribution  should  be  examined  to 
ensure  that  they  are  appropriate. 

A.  In  performing  these  steps  for  the  analysis  of  a  gravity  dam,  this 
phase  of  the  studv  is  limited  to  developing  a  method  to  analyze  the  deflec¬ 
tions  and  sL  reuses  of  the  gravity  concrete  structure  only.  Interaction  be¬ 
tween  the  structure  and  foundation  is  not  considered  at  this  time,  however  the 
reservoir  is  considered.  The  program  selected  in  Step  a  was  GTSTRUDL*  since 
it  is  supported  bv  t'e  vendor,  and  currently  widely  used  by  the  Corps.  Also, 
CToTRlTH.  is  ’-cpr-se-'tnt  ive  of  a  general-purpose  finite  element  program. 

PAST  :  :  •  ‘bis  report  presents  an  example  of  Steps  b-e  in  prepara¬ 

tion  for  the  anal  vs  is  <.f  a  gravity  dam.  An  actual  analysis  of  a  nonoverflow 
monolith  similar  to  the  Richard  B.  Russell  Dam  is  presented  in  PART  1IT.  The 
same  example  problem  is  then  analyzed  using  Chopra's  simplified  method  to 
illustrate  the  different  results  produced  by  each  method.  Part  IV  studies  the 
impose  of  including  a  foundation  with  regard  to  combined  static-dynamic 
stresses  within  a  gravity  dam. 


*  C7STRHDL  is  a  general-purpose  finite  element  program  owned  and  maintained 
by  the  GTTCES  Systems  Laboratory,  School  of  Civil  Engineering,  Georgia 
Institute  of  Technology.  Program  runs  used  in  this  report  were  made  on  the 
Control  Data  Corporation,  Cybernet  Computer  System. 
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PART  II:  FINITE  ELEMENT  DYNAMIC  ANALYSIS 
OF  A  SIMPLIFIED  STRUCTURE 


Selection  of  a  Simplified  Structure 

8.  The  simplified  cantilever  structure  representative  of  an  idealized 
gravity  dam  (Figure  I),  finite  element  models,  parameters,  along  with 


Figure  1.  Simplified  structure 


recommerdat ions  established  in  Phase  la  of  this  study  were  used  here.  Again, 
the  finite  element  runs  were  all  made  using  GTSTRUDL.  The  program  can  be  used 
in  the  analysis  of  the  static  and  dynamic  response  of  linear  two-  and  three- 
dimensional  (2-  and  3-D)  structural  systems.  The  element  used  was  the  "IPQQ" 
eight  node  isoparametric  quadratic  quadrilateral  element. 

Finite  Element  Models 


i in  it?  element  meshes 

-.  Three  different  models  previously  developed  in  the  Phase  la  report 
L.v  compare  the  convergence  characteristics  were  used  again.  The  various  mod¬ 
els  are  called  the  coarse,  fine,  and  very  fine  meshes  to  indicate  the  r2lative 
degree  of  refinement.  They  are  also  referred  to  as  Meshes  1,  2,  and  3  as 
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shown  below.  These  meshes  are  illustrated  in  Figures  2,  3,  and  4.  The  node 
and  elements  are  labeled  in  these  figures.  A  summary  of  the  meshes  is  pre¬ 
sented  below: 


Description 

Number  of  Nodes 

Number  of  Elements 
(All  IPQQ's? 

Mesh  1 

Coarse 

45 

10 

Mesh  2 

Fine 

149 

40 

Mesh  3 

Verv  fine 

537 

160 

Modeling  procedure 

10.  The  models  were  assumed  to  be  completely  restrained  along  their 
bases  and  to  be  in  a  state  of  plane  stress. 

Material  properties 

11.  The  weight  density  of  the  material  was  assumed  to  be  150  pcf.*  The 
modulus  of  elasticity  was  4,000,000  psi  with  a  Poisson  ratio  of  0.20. 

Dynamic  structural  properties 

12.  Dynamic  analysis  requires  the  same  input  to  describe  the  structural 
properties  as  does  a  static  analysis.  additional  requirements  are  that  the 
inertia  and  damping  of  the  structure  must  be  specified. 

13.  GTSTRUDL  will  automatically  compute  member/element  inertia  contri¬ 
butions  by  either  the  lumped  or  consistent  approaches.  The  member /element 
weight  densities  must  be  provided  via  the  CONSTANTS  command  or  the  MATERIALS 
command  prior  to  a  dynamic  analysis  if  automatic  computation  is  to  take  place. 
The  lumped  mass  approach  is  always  more  computationally  efficient  and  is  a 
reasonable  approximation  for  most  problems. 

14.  Damping  is  specified  in  GTSTRUDL  in  one  of  two  ways  depending  on 
whether  a  modal  superposition  or  direct  integration  transient  analysis  is  to 
be  performed.  In  this  study,  a  modal  analysis  will  be  performed,  thus  damping 
ratios  or  percent  damping  would  be  specified.  A  5  percent  damping  ratio  was 
assumed.  Had  the  stiffness  and  mass  matrices  been  input  via  the  MATRIX  com¬ 
mand,  damping  would  have  been  specified  by  proportional  damping  constants. 

15.  The  5  percent  damping  ratio  is  appropriate  for  a  mass  concrete  dam 
interacting  with  a  competent  rock  foundation  if  the  calculated  stress  levels 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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are  within  the  basic  strength  capacity  of  the  materials.  The  damping  ratio 
would  be  increased  If  the  stress  levels  go  above  the  basic  strength  capacity 
of  the  material. 

Dynamic  analysis 

16.  In  performing  a  dynamic  analysis,  the  GTSTRUDL  program  first  com¬ 
putes  the  mode  shapes  and  frequencies  of  the  structure.  These  results  are 
then  used  to  perform  a  time  history  analysis  or  a  response  spectrum  analysis, 
depending  on  the  choice  of  the  program  user. 

17.  The  response  spectrum  analysis  is  perhaps  the  most  common  technique 
used  in  many  design  offices.  The  computational  effort  required  for  a  time- 
history  analysis  is  often  prohibitive,  since  the  response  of  each  selected 
time  point  must  be  computed  and  stored  in  order  for  the  maximum  response  to  be 
identified.  Due  to  this  often  substantial  effort,  the  response  spectrum  anal¬ 
ysis  becomes  an  attractive  alternative  technique. 

18.  Response  spectrum  analysis  is  an  approximate  method  of  dynamic 
analysis.  It  uses  the  known  response  of  single  degree  of  freedom  systems  with 
the  same  natural  frequency  and  percents  of  critical  damping  as  the  modes  of 
vibration  of  the  structure  being  analyzed,  when  it  is  subjected  to  the  same 
transient  loading. 

19.  Once  the  maximum  response  of  each  mode  is  obtained,  the  maximum 
total  response  could  be  obtained  by  adding  the  maximum  response  of  each  mode 
since.  In  general,  however,  different  modes  will  attain  their  maximum  values 
at  different  times.  Therefore,  the  superposition  of  the  modal  maximums  will 
be  an  upper  bound  on  the  actual  total  response  and  will  significantly  over¬ 
estimate  the  response  for  many  cases. 

20.  GTSTRUDL  currently  computes  response  spectra  maximum  responses  by 
combining  the  modal  responses  by  seven  different  approaches.  ETL  1110-2-303 
(Department  of  the  Army  1985)  recommends  the  use  of  the  Complete  Quadratic 
Combination  (CQC)  method  (Der  Kiureghian  1980).  The  CQC  method  degenerates  to 
the  better  known  Square-Root-of-the-Sum-of-Squares  (SRSS)  method  for  simple, 

2-D  systems  in  which  the  frequencies  are  well  separated.  Combining  modal  max¬ 
ima  by  the  SRSS  method  can  dramatically  overestimate  or  significantly  under¬ 
estimate  the  dynamic  response,  especially  for  3-D  structures. 

Design  earthquake 

21.  Professor  H.  B.  Seed's  design  earthquake  response  spectrum  (Fig¬ 
ure  5)  (Seed,  Ugus,  and  Lysmer  1974),  scaled  to  0.25  g  peak  ground  acceleration 
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Figure  5.  Response  spectrum  (after  Seed,  Ugus, 
and  Lysmer  1974) 

with  5  percent  damping,  was  used  to  represent  the  earthquake  for  this  example 

22.  In  an  actual  analysis,  the  analyst  would  have  to  do  a  geological 
and  seismological  investigation  of  the  dam  site.  The  objective  of  the  inves¬ 
tigation  would  be  to  establish  the  controlling  earthquake,  and  the  correspond 
ing  ground  motions  to  be  used  in  the  study. 

Hydrodynamic  effects 

23.  It  has  long  been  understood  that  the  inertial  resistance  of  the 
water  in  a  reservoir  has  an  important  influence  on  the  earthquake  response  of 
concrete  dams.  This  study  considered  a  reservoir,  or  hydrodynamic  head  of 
170  ft  on  the  upstream  face  of  the  dam.  Except  in  the  case  of  the  EADHI 
(Chopra  and  Chakrabarti  1974)  and  EAGD-84  (Fenves  and  Chopra  1984)  codes, 
finite  element  models  usually  use  the  concept  of  an  "added  mass"  or  "virtual 
mass"  of  water  moving  with  the  dam  to  represent  the  hydrodynamic  interaction 
effect.  There  are  several  methods  of  approximating  this  effect.  The  one  cho 
sen  for  this  study  is  an  extension  of  the  Westergaard  method  that  was 
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originally  developed  for  gravity  dams  (Department  of  the  Army  1958, 

Westergaard  1933)  .  Westergaard  reasoned  that  the  added  mass  would  produce  the 
same  net  effect  on  lateral  loads  as  the  parabolic  hydrodynamic  pressure  dis¬ 
tribution.  The  effect  is  approximated  by  determining  and  attaching  added 
masses  to  the  fa^e  of  the  dam.  This  increased  mass  results  in  increased  iner¬ 
tial  resistance  to  the  motion  of  the  structure  when  an  earthquake  is  applied 
and  is  intended  to  simulate  the  actual  resistance  to  the  motion  of  the  struc¬ 
ture  caused  by  the  water  mass.  Calculation  of  the  added  mass  may  be  found  in 
Appendix  A. 

24.  Although  research  has  shown  the  Westergaard  added  mass  formulation 
provides  a  convenient,  simple  means  for  representing  reservoir  interaction  in 
the  dynamic  analysis  of  gravity  dams,  there  are  limitations  that  should  be 
noted.  The  underlying  assumptions  of  Westergaard ' s  work  are  that  the  dam  is 
rigid  and  the  water  is  incompressible.  Chopra  (1970)  has  shown  that  flexibil¬ 
ity  of  the  dam  and  compressibility  of  water  are  very  important  considerations 
in  the  dynamic  response,  hence  the  development  of  the  FAD^I  and  EAGD-84  code? 
(Chopra  1970).  The  main  drawback  to  these  programs  is  that  the  only  dynamic 
input  applicable  is  an  acceleration  time-history  record.  He  has  therefore 
developed  a  simplified  response  spectrum  analysis  procedure  (Chopra  1978)  for 
use  in  the  analysis  of  nonoverflow  concrete  gravity  dams  and  uncontrolled  (un¬ 
gated)  spillwTay  monoliths  which  can  be  modeled  for  2-D  analysis.  The  princi¬ 
ples  involved  are  those  basic  to  structural  analysis  by  response  spectrum 
methods.  The  method  represents  the  hydrodynamic  interaction  effects  by  an 
added  mass  of  water  which  moves  with  the  dam.  However,  unlike  Westergaard ' s 
added  mass,  this  mass  is  dependent  "on  the  frequency  and  shape  of  the  funda¬ 
mental  mode  of  vibration  of  the  dam,  and  the  effects  of  interaction  between 
the  flexible  dam  and  water,  considering  its  compressibility,  on  the  fundamen¬ 
tal  frequency  of  the  dam." 

Input 

25.  To  illustrate  the  use  of  GTSTRUDL  to  solve  a  problem  such  as  this 
simplified  structure,  the  input  file  used  for  the  Mesh  2  model  may  be  found  in 
Appendix  B. 
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Discussion  of  Results 


Mode  shapes  and  frequencies 

2 fi.  The  center-line  deflection  of  the  first  four  mode  shapes  of  the 
rectangular  beam  are  presented  in  Figure  6.  The  actual  deflected  shapes  are 
shown  in  Figures  7  to  10,  respectively.  GTSTRUDL  results  for  the  frequencies 
were  compared  to  Timoshenko  beam  theory  and  elementary  beam  theory  results 
(Table  2).  The  equations  used  to  generate  the  results  for  the  closed-form 
solutions  a-e  presented  in  Appendix  C.  Timoshenko  beam  theory  included  the 
effects  of  transverse  shear  and  rotary  inertia  which  is  significant  due  to  the 
relatively  short  and  deep  characteristics  of  the  structure. 


Figure  6.  Center-line  deflection  of  mode  shapes 
for  simplified  structure 
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Figure  9.  Third  mode  shape  of  simplified 
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Figure  10.  Fourth  mode  shape  of  simplified  structure 


Table  2 


Comparison  of  Frequencies  for  Simplified  Structure 


CTSTRUDL  (Mesh  2)  Includes  Added  Mass 


Fixed  Elementary  Timoshenko  (fixed  base)  GTSTRDDLt 


Mode 

Base 

Pin  and  Roller* 

Beam** 

Beam** 

Mesh  1 

Mesh  2 

Mesh  3 

1 

3.72 

3.63 

4.65 

4.11 

3.34 

3.34 

3.34 

2 

15.04 

13.52 

29.  15 

16.90 

12.90 

13.04 

13.09 

3+t 

15.37 

15.00 

— 

— 

15.03 

15.04 

15.04 

4 

32.67 

28.19 

81.59 

35.90 

25.82 

26.11 

26.25 

*  Although  a  fixed  base  should  be  used  in  an  analysis,  a  pin  and  roller  base 
was  also  run  for  comparison. 

**  Appendix  C  for  these  results. 

t  Frequencies  of  dam  including  effects  of  stored  water.  Other  frequencies 
do  not  include  any  stored  water  effects, 
i"*'  This  is  predominately  an  axial  mode,  Figure  6. 

27.  The  second  part  of  lable  2  indicates  the  deviation  in  frequencies 
due  to  the  various  model  refinements. 

Comparison  of  models 
for  the  dynamic  analysis 

28.  A  comparison  of  the  deflection  results  along  the  height  of  the 
rectangular  beam  is  presented  in  Table  3.  The  difference  in  the  maximum  dis¬ 
placement  is  approximately  3  percent  which  indicates  a  reasonable  agreement 
between  the  three  models. 

Table  3 

Comparison  of  Transverse  Deflections  Along  the  Simplified  Structure 


Distance  from 

Deflection,  in. 

Rase,  ft 

Mesh  1 

Mesh  2 

Mesh  3 

18.5 

0.021 

0.022 

0.022 

37.0 

0.064 

0.066 

0.066 

55.5 

0.122 

0.125 

0.  126 

74.0 

0.  193 

0.196 

0.197 

92.5 

0.273 

0.276 

0.277 

111.0 

0.358 

0.362 

0.363 

129.5 

0.446 

0.451 

0.452 

148.0 

0.536 

0.541 

0.543 

166.5 

0.625 

0.630 

0.632 

C*. 

LT\ 

QO 

0.711. 

('.716 

0.718 
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29.  A  comparison  of  the  SXY  (shear)  and  SYY  (normal)  stresses  was  made 
at  two  locations  along  the  height  of  the  simplified  structure.  The  results 
for  these  stresses  for  the  various  meshes  at  a  height  of  37.0  ft  above  the 
base  are  presented  in  Table  4,  w’hile  those  at  a  height  of  111.0  ft  are  in 
Table  5. 

30.  In  theory,  the  finer  the  mesh  the  more  correct  the  solution. 
Therefore,  Mesh  3  should  be  the  most  correct  solution.  The  results  indicate 
that  Meshes  1  and  2  are  converging  to  this  solution.  For  Mesh  1,  t!  differ¬ 
ence  from  Mesh  3  in  the  maximum  stresses  is  16  percent  for  shear  and  3  percent 
for  the  normal  stress.  Mesh  2  is  within  1  percent  in  both  cases  of  the  re¬ 
sults  of  Mesh  3.  Contour  plots  were  obtained  for  the  SXY  an  SYY  stress  compo¬ 
nents  for  Mesh  2  and  are  shown  in  Figures  11  and  12. 

31.  The  analyst,  in  addition  to  deciding  on  the  degree  of  mesh  refine¬ 
ment  so  as  to  achieve  sufficient  accuracy,  may  need  to  consider  cost.  The 
relative  costs  for  the  computer  runs  for  Meshes  1,  2,  and  3  was  $0.81,  $3.31, 
and  $23.25,  respectively.  (These  costs  are  from  runs  of  GTSTRUDL  on  the  Con¬ 
trol  Data  Cooperation  Cybernet  Computer  Service  and  should  be  used  only  as  a 
relative  measure.)  While  Mesh  3  should  produce  a  more  accurate  solution,  the 
additional  cost  does  not  appear  to  be  justified.  Mesh  2  provides  an  accept¬ 
able  solution,  balancing  both  cost  and  accuracy. 

Model  truncation  effects 

32.  An  investigation  was  made  to  determine  the  effects  of  using  various 
numbers  of  structural  vibration  modes.  Analyses  were  made  using  the  fundamen¬ 
tal  mode  and  subsequently  increasing  the  number  of  modes  until  the  final  solu¬ 
tion  showed  convergence.  Results  in  Table  6  show  no  difference  between  modes 
three  and  four,  thus  indicating  that  only  the  first  three  modes  are  necessary 
to  establish  complete  convergence.  Results  also  show  that  there  is  only  a 
minimal  difference  in  the  response  using  one  mode  and  the  computed  combined 
response  maximum  as  represented  by  the  four-mode  analysis.  This  indicates 
that  the  first,  or  fundamental,  mode  of  vibrations  participation  is  predomi¬ 
nate  in  obtaining  the  response-spectra  maximum  in  this  particular  dynamic 
analysis . 
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Table  4 

Comparison  of  Dynamic  Stresses  along  the 
Simplified  Structure  at  Height  of  37  Ft 


SYY,  psi 

SXY ,  psi 

X,  ft 

Mesh  1 

Mesh  2 

Mesh  3 

Mesh  1 

Mesh  2 

Mesh  3 

-40.0 

639 

612 

618 

20 

7 

5 

-30.0 

470 

468 

60 

62 

-20.0 

304 

313 

313 

97 

111 

107 

-10.0 

157 

156 

127 

129 

0.0 

8 

4 

4 

178 

141 

137 

10.0 

157 

157 

127 

129 

20.0 

305 

314 

314 

96 

111 

105 

30.0 

470 

468 

60 

62 

40.0 

638 

611 

618 

20 

5 

1 

Comparison  of 

Table  5 

Dynamic  Stresses  along 

the 

Simplified  Structure  at  Height  of  111 

Ft 

SYY,  psi 

SXY,  psi 

x,  ft 

Mesh  1 

Mesh  2 

Mesh  3  Mesh  1 

Mesh  2 

Mesh  3 

-40.0 

185 

184 

184  19 

8 

7 

o 

o 

m 

1 

145 

145 

39 

41 

-20.0 

102 

101 

101  65 

72 

69 

-10. 0 

52 

52 

83 

85 

0.0 

3 

3 

3  108 

93 

90 

10.0 

52 

52 

82 

83 

20.0 

102 

101 

} — ■ 

o 

Oh 

72 

70 

30.0 

147 

147 

37 

39 

40.0 

188 

187 

187  18 

4 

1 
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Table  6 

Effect  of  Number  of  Modes  Used  in  Analysis  on 
Dynamic  Stresses  (Mesh  2)  at  Height  of  37  Ft 


X,  ft 

SYY , 

psi 

SXY , 

psi 

1  Mode 

2  Modes 

3  Modes 

4  Modes 

1  Mode 

2  Modes 

3  Modes 

4  Modes 

-40.0 

611 

612 

612 

612 

5 

6 

7 

7 

-30.0 

470 

470 

470 

470 

58 

60 

60 

60 

-20.0 

313 

313 

313 

313 

107 

111 

111 

111 

-10.0 

157 

157 

157 

157 

123 

127 

127 

127 

0.0 

0 

4 

4 

4 

136 

141 

141 

141 

10.0 

157 

157 

157 

157 

122 

127 

127 

127 

20.0 

313 

314 

314 

314 

107 

111 

111 

111 

30.0 

470 

470 

470 

470 

58 

60 

60 

60 

40.0 

611 

611 

611 

611 

5 

5 

5 

5 
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PART  III:  GRAVITY  DAM  EXAMPLE  PROBLEM 


Description  of  the  Problem 


33.  An  actual  earthquake  analysis  using  both  the  dynamic  finite  element 
response  spectrum  method  as  outlined  in  Part  II  and  Chopra's  simplified  re¬ 
sponse  spectrum  method  is  presented  to  demonstrate  the  procedures  and  illus¬ 
trate  the  different  results  produced  by  each  method. 

34.  The  structure  to  be  analyzed  is  a  nonoverflow  monolith  similar  to 
those  of  the  Richard  B.  Russell  Dam.  The  dam  is  185  ft  high  with  a  reservoir 
depth  of  170  ft.  Seed's  design  response  spectrum,  scaled  to  a  peak  horizontal 
ground  acceleration  of  0.25  g  with  5  percent  damping  (Figure  5)  (Seed,  Ugus, 
and  Lysmer  1974),  will  be  used  in  both  analysis. 

35.  The  geometry  of  the  nonoverflow  section  is  defined  in  Figure  13. 


EL  495' 
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Finite  Element  Method  Analysis 

36.  A  listing  of  the  input  for  the  GTSTRUDL  finite  element  analysis  of 
the  gravity  dam  may  be  found  in  Appendix  D,  pages  D-2,  3,  and  4. 

37.  Previous  results  indicated  that  a  mesh  with  four  elements  across 
the  base,  Figure  14,  was  a  reasonable  compromise  between  accuracy  and  cost. 

The  mesh  contained  36  elements  and  135  nodes.  The  monolith  was  assumed  to  be 
completely  restrained  along  the  base. 

38.  The  structure  was  loaded  by  hydrostatic  and  hydrodynamic  loadings 
starting  at  170  ft  above  the  base  and  a  self-weight  of  the  concrete  of 

150  pcf.  The  hydrostatic  pressures  were  input  as  uniform  edge  loads  on  the 
upstream  elements.  The  hydrodynamic  effect  was  approximated  by  attaching 
Westergaard ' s  (1933)  "added  masses"  to  the  upstream  face  nodes,  Table  7. 
Analysis 

39.  The  analysis  is  performed  in  two  parts.  The  static  (stiffness) 
analysis  and  dynamic  analysis  are  performed  separately.  These  results  are 
then  combined  to  give  the  final  results.  The  static  analysis  consisted  of  two 
load  cases:  (1)  hydrostatic  pressure  on  the  upstream  face  of  the  dam,  and 

(2)  self-weight  (dead  load)  of  concrete. 

Results  of  analysis 

40.  Results  of  the  independent  load  cases  were  obtained.  It  should  be 
noted  that  elements  incident  on  a  common  node  will  have  different  stresses  at 
the  same  node.  This  is  due  to  the  fact  that  continuity  of  stresses  is  not  en¬ 
forced  or  required  for  the  finite  elements  in  GTSTRUDL,  as  is  true  in  all 
other  major  finite  element  programs.  To  obtain  a  more  useful  representation 
of  the  stresses,  one  can  use  the  CALCULATE  AVERAGE  command.  To  compute  the 
weighted  average,  GTSTRUDL  sums  the  stresses  for  all  elements  incident  on  a 
given  node,  and  then  divides  the  sum  by  the  number  of  elements  which  are  inci¬ 
dent  on  the  node. 

41.  Using  the  COMBINE  command,  it  is  then  possible  to  combine  the  in¬ 
dependent  loading  conditions  to  obtain  a  final  result.  In  this  example,  one 
would  add  the  two  static  loading  cases  (loads  1  and  2)  to  obtain  a  total 
static  loading  response  (loading  combination  5).  The  static  loading  condition 
than  can  be  combined  with  the  dynamic  loading. 

42.  It  is  first  necessary  to  operate  on  the  dynamic  loading  (load  3)  to 
transform  the  results  into  the  form  of  a  static  loading  condition.  The  result 
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GRAVITY  dam 


nodes 

Figure  H.  Mesh  for  dam  monolith 


Table  7 

Structure  Loading  for  Gravity  Dam  Example 


Hydrostatic 

Added 

Y 

Y 

y 

Pressure 

Mass 

Node 

Elevation 

ft 

ft 

ft 

psf 

slugs/ft 

127 

495.00 

185.00 

122 

113 

480.00 

170.00 

0.00 

2.50 

0 

55 

108 

5.00 

7.50 

312 

231 

09 

470.00 

160.00 

10.00 

14.25 

624 

463 

94 

18.50 

22.75 

1,154 

761 

85 

453.00 

143.00 

27.00 

33.07 

1,685 

1,136 

80 

39.13 

45.19 

2,442 

1,580 

71 

428.75 

118.75 

51.25 

57.32 

3,198 

1,811 

66 

63.38 

69.44 

3,955 

2,012 

57 

404.50 

94.50 

75.50 

81.57 

4,711 

2,199 

52 

87.63 

93.69 

5,468 

2,367 

43 

380.25 

70.25 

99.75 

105.82 

6,224 

2,527 

38 

111.88 

117.94 

6,981 

2,675 

29 

356.00 

46.00 

124.00 

129.75 

7,738 

2,742 

24 

135.50 

141.25 

8,455 

2,793 

15 

333.00 

23.00 

147.00 

152.75 

9,173 

2,909 

10 

158.50 

164.25 

9,890 

3,021 

1 

310.00 

0.00 

170.00 

10,608 

1,551 

Notes : 

Y 

is  measured 

from  base 

of  dam  (el 

310) 

Y 

is  depth  below  water 

surface  (el 

480) 

•  * 

v  Is  depth  below  water  surface  to  midpoint  between  nodes. 
Ce  and  Mi  are  defined  in  Appendix  A. 

Ce  =  - -  51-—  ■—  =  51.54 

•f-  °-72  (uoo67°  i.o)2 


* 


_  2  «  51.54  x  170 
i  3  x  32.2 

All  elevations  (el) 
Vertical  Datum  (NGVD) 


=  13.91 


added  mass 


cited  herein  are  in  feet  referred  to  the  National  Geodetic 
of  1929. 
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is  called  a  pseudostatic  loading.  The  command  performs  this  function  by 
copying  the  results  of  the  specified  modal  combination  (CQC)  of  the  response 
spectrum  loading  condition  into  a  static  loading  condition  (loading  4).  This 
can  then  be  added  or  subtracted  with  the  static  loading  condition  to  obtain 
the  maximum  or  minimum  stresses  in  the  dam. 

43.  A  listing  of  the  output  can  be  found  in  Appendix  D.  Page  D-5  lists 
the  dams  first  four  frequencies.  Pages  D-6  through  D-21  give  the  stresses  for 
the  loading  conditions;  load  4  -  pseudostatic  loading,  load  5  -  static  load¬ 
ing,  and  load  6  -  static  plus  dynamic  loading.  Extracted  from  the  output  and 
presented  in  Table  8  are  the  vertical  (SYY)  stresses  for  the  monoliths 

Table  8 

Stresses  in  Dam 


Node 

SYY,  psi 

Static 

Dynamic 

Static  plus 
Dynamic 

1 

53 

313 

366 

9 

-41 

44 

3 

15 

-23 

275 

253 

23 

-60 

88 

29 

29 

-59 

260 

202 

37 

-89 

167 

79 

43 

-56 

261 

205 

51 

-78 

192 

114 

57 

-51 

260 

208 

65 

-56 

197 

141 

71 

-49 

257 

208 

79 

-33 

199 

166 

85 

-49 

273 

225 

93 

-13 

194 

181 

99 

-28 

188 

160 

107 

-29 

285 

256 

113 

-15 

79 

64 

121 

-14 

54 

40 

Note : 

Positive  -  Tension 

Negative  =  Compression 
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upstream  and  downstream  face  nodes.  Column  2  presents  the  weighted  average 
stress  at  '•F‘=i  various  "<->Hes  due  to  the  static  Tending.  In  Column  3,  the 
dynamic  stresses  are  given.  The  final  column  presents  the  combined  stresses, 
the  addition  of  the  static  and  dynamic  stresses.  This  would  represent  the 
maximum  tension  in  the  monolith  due  to  the  prescribed  earthquake  with  an  up¬ 
stream  pool  at  el  480  ft. 

44.  The  upstream  and  downstream  face  node  SYY  stress  component  along 
with  the  contour  plots  for  the  various  loading  conditions  are  shown  in  Fig¬ 
ures  15,  16,  and  17. 

45.  As  a  sidelight  to  the  above  example  problem,  one  additional  analy¬ 
sis  was  made.  As  previously  mentioned,  the  inertial  resistance  of  the  water 
in  a  reservoir  has  an  influence  on  the  earthquake  response  of  a  dam.  But,  how 
great  is  the  influence?  To  answer  this,  an  additional  computer  run  was  made 
without  the  hydrodynamic  effects  (attached  "added  masses").  Figure  18  shows 
the  resulting  face  node  SYY  stress  components.  Comparing  this  with  Figure  16 
one  can  see  the  higher  stresses  in  the  dam  due  to  dam/reservoir  interaction. 

Chopra  Simplified  Analysis 

46.  The  more  rigorous  dynamic  analysis  of  gravity  dams  is  obtained  by 
using  finite  element  computer  programs.  Due  to  the  capability  of  these 
programs  to  model  the  horizontal  and  vertical  structural  deformations  of  the 
dam,  to  model  the  exterior  and  interior  concrete,  and  to  include  the  response 
of  the  higher  modes  of  vibration,  the  interaction  effect  of  the  foundation  and 
any  surrounding  soil,  and  the  horizontal  and  vertical  components  of  the  ground 
motion,  some  amount  of  specialized  training  is  required  to  use  them  effec¬ 
tively.  Chopra's  simplified  analysis  procedure  is  a  compromise  on  that 
complexity . 

47.  Using  a  set  of  standard  curves  for  the  fundamental  mode  shape,  the 
ratio  of  fundamental  period  of  the  dam  with  and  without  stored  water,  and  the 
variation  of  hydrodynamic  pressure  over  the  depth  of  the  water,  one  can  calcu¬ 
late  a  set  of  equivalent  static  lateral  loads.  These  forces  are  considered  to 
act  separately  in  the  upstream  and  downstream  directions,  and  their  effects 
added  to  the  effects  of  all  other  design  loads. 
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48. 


The  computation  of  the  earthquake  forces  is  carried  out  as  follows: 

a.  Compute  Tg  ,  the  fundamental  natural  period  of  vibration  of 

the  dam,  in  seconds,  without  the  influence  of  stored  water, 
from 


T 


S 


1.4  x  H 
_ s 

/e 


CD 


where 

H  =  height  of  dam,  feet 

s 

E  =  modulus  of  elasticity  of  concrete,  psi 

b.  Compute  ,  the  fundamental  natural  period  of  vibration  of 

the  dam,  in  seconds,  with  stored  water 


T  =  R  T  (2) 

s  Is 

where  R  =  period  ratio  determined  from  Figure  19. 

c.  Compute  ,  the  ratio  of  the  fundamental  resonant  period  for 

the  impulsive  hydrodynamic  pressure  to  T  from 

.  (3) 

T 

s 


where 

C  =  velocity  of  sound  in  water  =  4,720  ft/sec 
H  =  depth  of  water  (feet)  in  Figure  19 

d.  Compute  f^Cy)  *  the  lateral  earthquake  forces  over  the  height 
of  the  dam,  including  hydrodynamic  effects,  from 


A,  x  S  (T  )  _ 

fg(y)  =  - - - Jws(y)<Ky)  +  gPjCy)]  CM 


where 


=  scaling  constant,  with  assumed  value  of  4 
S  (T  )  =  a  spectral  acceleration  at  period  of  vibration 

3.  S 

Tg  ,  from  the  design  response  spectrum 
g  =  acceleration  due  to  gravity 
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<Ky)  =  fundamental  mode  shape  factor  from  Figure  20 

gp^(y)  “  pressure  distribution  factor  from  Figure  21  cor¬ 
responding  to  R2  and  multiplied  by  the  quantity 

(H/H  )2 
s 

Compute  f(y)  ,  the  lateral  earthquake  forces  without  hydro- 
dynamic  effects  from: 

S  (T  ) 

f(y)  =  wg(y)<Ky)  (5) 

where 


S  (T  )  =  spectral  acceleration  at  period  of  vibration  T 
Si  s  s 


<a  cc 
b-  uj 


I  § 


Figure  19.  Illustration  of  how  to  determine  standard 
values  for  ,  the  ratio  of  fundamental  vibration 

periods  of  the  dam  with  and  without  water,  for  a 
given  H/H  ratio  and  modulus  of  elasticity  (after 
S  Chopra,  1978) 


MOOE  SHAPE.  0 


Figure  20.  Illustration  of  how  to 
determine  the  standard  fundamental 
period  and  mode  shape  for  a  given 
y/H^  ratio  (after  Chopra  1978) 


Figure  21.  Illustration  of 
how  to  determine  th£  stan¬ 
dard  variation  of  gp^  over 

the  depth  of  water  for  H/H 

=  1  and  various  values  of  S 

R  where  W  is  the  unit 
z 

weight  of  water  (after 
Chopra  1978) 
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49.  The  equations  for  distributed  lateral  force  (Equations  4  and  5), 
can  be  integrated  between  appropriate  limits  to  yield  concentrated  forces 
which  are  then  applied  statically  to  the  dam.  Stress  computations  are  then 
carried  out  as  with  other  design  loads. 

Computation  of  earthquake  forces 

50.  For  this  analysis,  the  nonoverflow  monolith  of  the  gravity  dam  is 
divided  into  nine  sections  (or  horizontal  slices).  The  elevation  of  each  sec¬ 
tion  is  equal  to  the  corner  node  elevations  of  the  elements  used  in  the  pre¬ 
vious  FEM  analysis.  Thus,  section  1  is  at  an  elevation  equal  to  that  for 
nodes  113  through  121  (see  Figure  14),  section  2  is  at  nodes  99  through  107, 
section  3  at  nodes  85  through  93,  etc. 

51.  Steps  in  the  computation  of  the  earthquake  forces  are  as  follows: 

a.  For  E  =  4.0  *  10^  psi  and  dam  height  H  =  185  ft,  from 
Equation  1, 


T 

s 


1.4  x  185 
4.0  x  106 


0.13  sec 


b.  From  Figure  18, 

R  =  1.23  for  E  =  4.0  x  io6  psi 

and 


H 

H 


170 

185 


=  0.92 


From  Equation  2, 


T  =  1.23(0.13)  =  0.16  sec 
s 


£.  From  Equation  3, 


_  4.0  x  (170/4,720) 
2  "  0.16 


0.90 


d.  Equation  4  was  evaluated  along  each  level  throughout  the  height 

Sa(TsVg]  -  0. 


of  the  dam,  by  substituting 


A1  =  4 


and 


.Hi. 

[' 


.63 


(from  Figure  5)  for  Tg  =  0.16  sec,  by  computing  the  weight  of 
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the  dam  per  unit  height,  wg(y)  »  from  the  section  dimensions 

(Figure  13)  and  the  unit  weight  of  concrete,  and  by  substitut¬ 
ing  for  <J>  (y)  from  Figure  19  and  gPj(y)  ,  for  Figure  21. 

The  distributed  lateral  (earthquake)  forces  f  (y)  and  equiv¬ 
alent  static  load  are  listed  in  Table  9  and  pictorially  shown 
in  Figures  22  and  23. 

Computation  of  stresses 

52.  At  this  point,  simple  beam  theory 


could  be  used  to  calculate  the  stresses  at  each  horizontal  section.  As  was 
done  with  the  earthquake  forces,  the  distributed  gravity  and  hydrostatic 
forces  are  replaced  by  concentrated  loads.  The  direct  and  bending  stresses 
are  then  computed  based  on  section  properties  (see  Table  10)  at  each  horizon¬ 
tal  slice. 

53.  Tables  11  and  12  show  the  results  of  the  stress  calculations.  The 
earthquake  loads  have  been  applied  both  upstream  and  downstream  and  the  re¬ 
sults  combined  with  the  static  loads.  Figure  24  shows  the  static  stresses 
(SYY)  while  Figure  25  shows  the  maximum  tensile  stresses  (SYY),  a  result  of 
combining  the  static  and  dynamic  stresses. 

54.  An  alternative  to  using  simple  beam  theory  to  calculate  stresses 
would  be  to  perform  a  static  finite  element  analysis.  This  approach  has  been 
implemented  by  Cole  and  Cheek  (Technical  Report  SL-86-44,  Department  of  the 
Army)  in  a  new  user-friendly  computer  program.  The  program  was  developed 
using  the  finite  element  methods  of  analysis  to  determine  the  dam's  inertial 
response  along  with  Chopra's  simplified  procedure  for  estimating  the  hydro- 
dynamic  loading.  The  program  is  menu  driven,  allowing  for  ease  of  use  by  the 
novice.  The  only  input  required  is  the  dam's  geometry  and  appropriate  re¬ 
sponse  spectrum.  Output  consists  of  principal  surface  stress  values  at 
selected  elevations  for  both  the  upstream  and  downstream  faces  of  the  dam.  A 
run  using  this  program  is  included  in  Appendix  E.  The  results  are  tabulated 
in  Appendix  E,  page  E4  and  shown  in  Figure  26. 

Comparison  of  Procedures 

55.  Corresponding  calculated  stresses  are  plotted  for  each  method  of 
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SCALE 

L-  ■  .1  .~-.J _ _  ■  \  . 

2  1  Ox  104  LB/FT 

Figure  22.  Dynamic  force  (inertia  +  water)  on  gravity  dam 
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3  2  1  0  x  10s  LBS 

Figure  23.  Dynamic  load  (inertia  +  water)  on  gravity  dam 


Table  10 

Section  Properties  of  Gravity  Dam 


Base 

Base 

Width 

W 

A  2 

s  , 

Section 

Elevation 

ft 

lb 

in. 

in. 

1 

480.00 

17.00 

38,250 

2,448 

83,232 

2 

470.00 

17.00 

63,750 

2,448 

83,232 

3 

453.00 

28.33 

121 ,546 

4,080 

231,254 

4 

428.75 

46.52 

257,679 

6,699 

623,264 

5 

404.50 

64.71 

465,435 

9,318 

1,205,967 

6 

380.25 

82.89 

733,883 

11,936 

1.978,777 

7 

356.00 

101.08 

1,068,478 

14,556 

2,942,544 

8 

333.00 

122.52 

1,454,188 

17,643 

4,323,211 

9 

310.00 

143.25 

1,912,641 

20,628 

5,909,922 

analysis;  finite  element  analysis  using  Westergaard ' s  "added  mass"  (FEM) , 
Chopra's  simplified  method  using  simple  beam  theory  (CSM/SBT) ,  and  Chopra's 
simplified  method  using  finite  element  analysis  (CSM/FEM) .  Comparisons  are 
presented  for  both  the  upstream  (Figure  27)  and  the  downstream  (Figure  28) 
faces . 
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Table  11 

Vertical  Stresses  in  Gravity  Dam  with  Earthquake  In 
Downstream  Direction  Using  Simplified  Method 

Static  Stresses 
+ 

Static  Stresses  Dynamic  Stresses  Dynamic  Stresses 

psi  psi  psi 


Sect! on 

Flevation 

UDStream 

Downstream 

Upstream 

Downstream 

Upstream 

Downstream 

1 

480.00 

-16 

-16 

93 

-93 

77 

-109 

2 

470.00 

-24 

-28 

257 

-257 

233 

-285 

3 

453.00 

-6 

-54 

275 

-275 

269 

-329 

4 

428.75 

-30 

-46 

283 

-283 

253 

-329 

5 

404.50 

— 

-56 

299 

-299 

255 

-355 

6 

380.25 

-50 

-72 

311 

-311 

261 

-383 

7 

356.00 

-52 

-94 

318 

-318 

266 

-412 

8 

333.00 

-62 

-102 

299 

-299 

237 

-401 

9 

310.00 

-70 

-116 

287 

-287 

217 

-403 

Note:  Positive  =  Tension. 

Negative  =  Compression. 


Table  12 

Vertical  Stresses  in  Gravity  Dam  with  Earthquakes  in 
Upstream  Direction  Using  Simplified  Method 


Static  Stresses 
+ 


Static 

Stresses 

Dynamic 

Stresses 

Dynamic 

Stresses 

psi 

psi 

psi 

Section 

Elevation 

Upstream 

Downstream 

Upstream 

Downstream 

Upstream 

Downstream 

1 

480.00 

-16 

-16 

-  93 

93 

-109 

77 

2 

470.00 

-24 

-28 

-257 

257 

-281 

229 

3 

453.00 

-6 

-54 

-275 

275 

-281 

221 

4 

428.75 

-30 

-46 

-283 

283 

-313 

237 

5 

404.50 

-44 

-56 

-299 

299 

-343 

243 

6 

380.25 

-50 

-72 

-311 

311 

-261 

239 

7 

356.00 

-52 

-94 

-318 

318 

-370 

224 

8 

333.00 

-62 

-102 

-299 

299 

-361 

197 

9 

310.00 

-70 

-116 

-287 

287 

-357 

171 

Note : 

Positive  ■ 

Tension. 

Negative  »  Compression. 
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Figure  24 .  Static  stresses  (SYY)  in  gravity  dam  from  simplified  method 
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Figure  25.  Maximum  tensile  stresses  (SYY)  in  gravity  dam  from 

simplified  method 
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PART  IV:  FOUNDATION  EFFECTS 


Finite  Element  Model 


56.  The  purpose  of  this  portion  of  the  study  in  to  determine  the  neces¬ 
sity  of  including  a  foundation  in  the  finite  element  model  used  for  the  dy¬ 
namic  analysis  of  a  gravity  dam.  The  same  gravity  dam  geometry  and  mesh  size 
as  described  in  the  Part  III  FEM  analysis  were  used  for  this  study  with  the 
addition  of  a  foundations  block  (Figure  29).  Limits  on  the  width  and  depth  of 
a  foundation  block  were  evaluated  in  the  Phase  lb  report  (in  preparation). 

The  model  used  in  this  study  included  foundation  material  in  both  the 
upstream  and  downstream  directions  equal  to  the  base  width  of  the  dam 
and  a  depth  equal  to  1.5  times  the  base  width.  The  resulting  mesh  contains 
68  elements  and  247  nodes. 


Y 
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57.  The  foundation  boundary  was  assumed  fixed  at  the  corners  ("nodes  1 
and  17)  and  free  to  slide  along  its  base  and  sides.  Negligible,  mass  was  given 
to  the  foundation  to  avoid  wave  propagation  effects  at  this  level. 


Stress  Comparison 


58.  In  order  to  assess  the  effect  of  the  foundation,  the  dam  was  ana¬ 
lyzed  five  times  varying  only  the  modulus  of  elasticity  of  the  foundation 
material.  The  various  Er/Ec  ratios  (the  ratio  of  the  modules  of  elasticity  of 
the  foundation  (Er)  to  the  dam  (Ec))  were:  0.05,  0.25,  1.00,  1.75,  and  3.00. 
An  analysis  assuming  an  Er/Ec  ratio  of  infinite  (°°)  is  equivalent  to  the 
analysis,  as  performed  in  Part  III,  assuming  the  dam  is  completely  restrained 
along  its  base. 

59.  Comparison  of  the  first  mode  frequency  and  period  for  the  various 
Er/Ec  ratios  is  presented  below: 


Er/Ec 

Frequency 

(cyc/sec) 

Period 

(sec/cy 

Infinite 

6.042 

0.165 

3.00 

5.543 

0.180 

1.75 

5.256 

0.190 

1.00 

4.836 

0.207 

0.25 

3.353 

0.298 

0.05 

1.731 

0.578 

The  addition  of  a  foundation  to  a  model,  regardless  of  its  size,  will  decrease 
its  frequency  and  increase  its  period.  This  substantiates  that  in  addition, 
as  the  foundation's  modules  of  elasticity  decrease,  it  will  further  reduce  the 
frequency  and  increase  the  period. 

60.  The  upstream  and  downstream  dam  face  node  SYY  stress  components 
along  with  the  contour  plots  for  the  static,  dynamic,  and  combined  static  and 
dynamic  analysis  are  obtained  for  the  various  Er/Ec  ratios  and  .  lO’-m  in  Fig¬ 
ures  30  through  44.  Previously  presented  in  Figures  15,  16,  and  17  were  the 
stress  components  and  'ontour  plots  for  the  analysis  assuming  a  completely 
restrained  base  (Er/Ec  =  infinite).  Combined  results  for  the  SYY  stress  com¬ 
ponents  are  presented  in  Figures  45  through  50. 

61.  Figures  45  through  46  show  that  as  the  foundation's  modulus  of 
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DEPTH  DEPTH 


DEPTH 


Figure  35.  Maximum  tensile  (static  plus  dynamic)  stresses  (SYY) 

for  Er/Ec  ratio  0.25 


ELEVATION 


elevation  elevation 


-22  -66  3  20 

Figure  38.  Maximum  tensile  (static  plus  dynamic)  stresses  (SYY) 

for  Er/Ec  ratio  of  1.00 


-96  -105  -'08 


Figure  39.  Static  stresses  (SYY)  for  Er/Ec  ratio  of  1.75 


DEPTH 


Figure  46.  Downstream  face:  Static 
stresses  (SYY) 
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STRESSES  -  PS1 

Figure  48.  Downstream  face:  Dynamic  tensile  stresses  (SYY) 
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Figure  49.  Upstream  face:  Static  +  dynamic  tensile  stresses  (SYY) 


0  SO  100  ISO  700  750  300  350 

STRESSES  PSI 


Figure  50.  Downstream  face 


Static  +  dynamic  tensile  stresses  (SYY) 
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elasticity  decreases,  the  static  stresses  (SYY)  below  midheight  of  the  dam 
increase.  This  was  confirmed  in  the  Phase  lb  report  (in  preparation). 

62.  Figures  47  and  48,  showing  the  dynamic  stresses,  are  net  quite 
clear  cut.  To  understand  the  results,  one  must  first  look  at  the  response 
spectrum.  Principally,  the  first  mode  of  vibration  results  in  the  largest 
participation  in  the  response  spectra  maximum.  Therefore,  the  model  whose 
fundamental  mode  shapes  period  results  in  the  largest  corresponding  maximum 
response  values  will  generally  have  the  highest  dynamic  stresses. 

63.  Seed's  smoothed  response  spectrum  (Figure  5)  was  used  in  this  anal¬ 
ysis.  Figure  51  shows  the  actual  digitized  response  spectrum  used  for  the 
computer  programs  input.  The  first  mode  periods  associated  with  the  various 
Er/Ec  ratios  are  plotted  on  this  response  curve  along  with  the  corresponding 
maximum  response  values.  It  so  happens,  in  this  example,  that  the  model  with 
an  Er/Ec  of  infinite  results  in  the  peak  maximum  response  of  0.63g.  As  the 
Er/Ec  ratio  decreases,  so  does  the  response  acceleration.  This  would  indicate 


Figure  51.  Digitized  response  spectrum 
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that  the  analysis  using  the  Er/Ec  ratio  of  infinite  would  result  in  the  high¬ 
est  dynamic  stresses,  and  stresses  would  decrease  as  the  Er/Ec  decreased. 

64.  In  this  particular  case,  results  shown  in  Figures  47  and  48  do  not 
quite  bear  this  out.  The  maximum  stresses  resulted  from  an  Er/Ec  ratio  of 
1.00.  Er/Ec  ratios  of  1.75,  3.00,  and  infinite  resulted  in  slightly  decreas¬ 
ing  stresses.  This  is  opposite  of  what  would  be  expected  and  must  be  attrib¬ 
uted  to  foundation-structure  interaction.  Stresses  for  the  four  highest  Er/Ec 
ratios  were,  however,  generally  within  5  percent  of  each  other. 

65.  Er/Ec  ratios  of  0.25  and  0.05  did  result  in  lower  stresses,  as 
expected,  except  in  the  lower  portion  of  the  dam. 

66.  Figures  49  and  50  show  the  combined  static  and  dynamic  stresses. 
Comparison  of  the  various  plots  show  that  combined  stresses  for  Er/Ec  ratios 
above  0.25  are  within  25  percent  of  each  other  except  at  the  toe  of  the  dam. 
This  indicates,  for  this  particular  case,  the  finite  element  grid  need  not 
include  a  foundation  until  the  foundation  materials  modulus  of  elasticity  is 
less  than  25  percent  of  that  in  the  dam. 
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PART  V:  SUMMARY  AND  RECOMMENDATIONS 


Summarv 


67.  The  primary  objective  of  this  study  was  to  illustrate  an  approach 
for  performing  a  finite  element  response-spectrum  dynamic  analysis  of  a  grav¬ 
ity  dam.  Ihe  illustrations  should  serve  the  beginning  finite  element  analyst 
with  a  better  understanding  of  the  behavior  of  a  gravity  dam  which  is  subject 
to  dynamic  loading. 

68.  The  analysis  of  a  monolith  similar  to  those  of  the  Richard  B. 
Russell  Dam  determined  that  a  mesh  with  four  elements  across  the  base  was  a 
reasonable  compromise  between  accuracy  and  cost  in  the  dynamic  analysis.  The 
finite  element  grid  need  not  include  a  foundation  bTock  as  long  as  the  founda¬ 
tion  materials  modulus  of  elasticity  is  at  least  26  percent  of  that  in  the 
dam. 

69.  Results  of  the  comparison  between  a  FEM  analysis  using 
Westergaard 1 s  (1933)  "added  mass"  and  Chopra's  (1978)  "simplified  response 
spectrum"  method  showed  the  simplified  method  to  be  conservative,  but  not 
excessively  oO.  Results  indicate  the  simplified  method  could  be  used  to  make 
a  first-cut  estimate  of  the  surface  stresses. 

Recommendations 


70.  Conclusions  reached  in  this  report  were  based  on  studies  using  a 
single  monolith  size.  Varying  the  overall  dimensions  of  the  dams  may  signifi¬ 
cantly  alter  the  results.  This  report  should  illustrate  to  the  engineer  the 
importance  of  making  verification  studies  to  ensure  the  use  of  a  proper  mesh 
size  and  the  necessity  of  including  a  foundation  from  which  usable  results  can 
be  obtained.  The  analyst  should  develop  and  analyze  finite  element  models  of 
simplified  structures,  then  extrapolate  the  information  gained  to  the  modeling 
of  the  real  structure. 
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APPENDIX  A:  MODELING  OF  HYDRODYNAMIC  EFFECTS 


1.  The  "added  mass"  applied  to  a  structure  to  simulate  the  hydrodynamic 
effects  can  be  computed  using  Westergaard ' s  formula,  EM  1110-2-2200  (Depart¬ 
ment  of  the  Army  1958),*  which  gives  the  hydrodynamic  pressure  at  a  depth  y 
below  the  water  surface  as: 

p  =  C^a/hy  lb/ft^  (1) 

P  =  f  Ce  I  (?) 

where 

p  =  hydrodynamic  pressure  at  depth  y  below  water  surface,  pounds  per 
square  foot 

P  =  total  pressure  to  depth  y  from  surface  using  the  parabolic, 
pounds  per  square  foot  approximation 

h  =  total  depth  of  water,  feet 

a  =  ratio  of  earthquake  acceleration,  a  to  g 

a  =  acceleration  due  to  the  earthquake,  feet  per  second  squared 

2 

g  =  gravitation  acceleration,  32.2  ft/sec 

Cg  =  a  factor  depending  principally  on  height  of  dam  and  the  earthquake 
vibration  period,  tg  ,  sec 

2.  Westergaard ' s  approximate  equation  for  Cg  ,  which  is  sufficiently 
accurate  for  all  usual  conditions,  in  pounds  per  cubic  foot  is: 

51  3 

Ce  =  — - -  ■■■  ib/ff5  (3) 

yj 1  '  0,72  (t^ot) 

Period  of  vibration,  t^  ,  is  usually  assumed  as  1  sec.  The  mass  per  unit 
area  to  be  added  to  the  face  of  the  dam  is  then  calculated  by  dividing  the 
pressure  by  the  acceleration,  a  .  This  gives: 

m  =  —  /hy  lb-sec  /ft  (4) 

*  References  cited  in  this  appendix  are  included  in  the  references  at  the  end 
of  the  main  text. 
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lb-sec^/ f 


where 

m  =  mass  per  unit  area  to  be  added  to  the  face  of  the  dam 
M  =  total  mass  to  depth  y 

The  total  mass  to  be  added  to  a  particular  area  on  the  face  of  the  dam  is  then 
found  by  integrating  this  quantity  over  the  area  under  consideration.  These 
added  masses  are  lumped  at  the  node  points  of  the  finite  element  grid  on  the 
face  of  the  dam.  This  gives: 


M. 

1 


2  Ce  0 . 5 
w  —  h 

3  g 


slugs/f  t 


(6) 


where 


M. 

x 


y 


2 


added  mass  to  be  applied  at  node 


i 


1 


pool  depth  to  the  midpoint  between  node  i  and  the  node  directly 
below 


y 


l 


pool  depth  to  the  midpoint  between  node  i  and  the  node  directly 
above 
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APPENDIX  B:  INPUT  FILE  FOR  MESH  2  MODEL 


STRUDL  'MESH  2'  'FINE  MESH' 

$  4  BY  10  MESH 
UNITS  KIPS  FEET 

GEN  9  JOI  ID  1  1  X  0.0  10.0  Y  0.0 
MOD  10  ID  14  Y  18.5 

GEN  5  JOI  ID  10  1  X  0.0  20.0  Y  9.25 
MOD  9  ID  14  Y  18.5 
TYPE  PLANE  STRESS 


GEN 

10 

ELE 

ID 

1 

4 

FROM 

1, 

14 

TO 

3 

TO 

17 

TO 

15 

TO 

2 

TO 

11 

TO 

16 

TO 

10 

GEN 

10 

ELE 

ID 

2 

4 

FROM 

3, 

14 

TO 

5 

TO 

19 

TO 

17 

TO 

4 

TO 

12 

TO 

18 

TO 

11 

GEN 

10 

ELE 

ID 

3 

4 

FROM 

5, 

14 

TO 

7 

TO 

21 

TO 

19 

TO 

6 

TO 

13 

TO 

20 

TO 

12 

GEN 

10 

ELE 

ID 

4 

4 

FROM 

7, 

14 

TO 

9 

TO 

23 

TO 

21 

TO 

8 

TO 

14 

TO 

22 

TO 

13 

STAT  SUPPORT  1  TO  9 
CONSTANTS 
E  576000.  ALL 
POISSON  0.20  ALL 
DEN  0.150  ALL 
ELEM  PROP 

1  TO  40  TYPE  ' IPQQ '  THICK  1.0 
DAMPING  0.05  4 
* 

UNITS  KIPS  FEET  SECONDS  CYCLES 

STORE  RESPONSE  SPECTRA  ACCELERATION  LIN  OS  PERIOD  LIN  'SEED'  DUMP 
$  ACCELERATION  ( FT/SEC AA2 )  US  PERIOD  (SEC) 


DAMPING  0.05 

FACTOR 

0.25 

31 .78 

.0050 

33.42 

.0260 

51.23 

.  0745 

59.70 

.0805 

73.99 

.  1260 

81.69 

.  1635 

72.61 

.2680 

68.68 

.2890 

57.06 

.3755 

53.19 

.4005 

40.22 

.5055 

37.22 

.  5415 

27.05 

.6980 

23.92 

.7955 

17.87 

1.0905 

16.00 

1 .1945 

9.98 

1 .6610 

8.82 

1.8030 

END  OF  RESPONSE  SPECTRUM 
$ 

UNITS  KIPS  FEET  SECONDS 
INERTIA  OF  JOI  LUMPED 
INERTIA  OF  JOINTS  MASS 
$  HYDRODYNAMIC  'ADDED  MASS' 


1 

TRANS 

X 

1.41 

10 

TRANS 

X 

2.76 

15 

TRANS 

X 

2.68 

2  4 

TRANS 

X 

2.60 

29 

TRANS 

X 

2.51 

B  3 


29  TRANS  X  2.51 
38  TRANS  X  2.42 
43  TRANS  X  2.33 
52  TRANS  X  2.23 
57  TRANS  X  2.13 
66  TRANS  X  2.03 
71  TRANS  X  1.92 
80  TRANS  X  1.80 
85  TRANS  X  1.67 
94  TRANS  X  1.54 
99  TRANS  X  1.39 
108  TRANS  X  1.22 
113  TRANS  X  1.02 
122  TRANS  X  0.77 
127  TRANS  X  0.36 
* 

UNITS  LBS  INCHES 

DYN  LOAD  1  'SEED  RESPONSE  SPECTRUM' 

SUPPORT  ACC 

TRANSLATION  X  FILE  'SEED' 

END  °F  DYN  LOAD 
$ 

$ 

E IGENPROBLEM  PARAMETERS 
SOLVE  USING  SUBSPACE  ITERATION 
NUMBER  OF  MODES  4 

PERFORM  NO  STRUM  SEQUENCE  CHECK  ORTHOGONALITY  CHECK 
TOLERANCE  E IGENVAL  l.E-4 
END 
$ 

DYN  ANAL  MODAL 

COMPUTE  DYNAMIC  DISPLACEMENTS  FORCES  STRESSES  MODAL  COMB  ALL 

PRINT  DYN  DATA 

OUTPUT  DECIMAL  4 

OUTPUT  BY  LOADING 

OUTPUT  FIELD  E 

LIST  DYNAMIC  EIGENVE 

CREATE  PSEUDO  STATIC  LOADING  2  'CQC  OF  LOADING  1'  AS  CQC  OF  LOADING  1 

DELETIONS  ;  LOAD  1 

LIST  DISPL  STRESSES 

CALCULATE  AVERAGE  STRESSES 

SAVE  DIRECT  'MESSAV1' 

FINISH 
END  OF  FILE 
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APPENDIX  C:  NATURAL  FREQUENCY  CALCULATIONS  OF  A  CANTILEVER  BEAM 


Elementary  Beam  Theory 


1.  This  discussion  concerns  natural  frequency  of  a  cantilever  beam 
based  only  on  the  elementary  engineering  theory  of  beam  bending  with  no  sec¬ 
ondary  effects  (Warburton  1964):* 


0.  132L 


3,516/  El\  1/2 

22.03  / E I \  ’/2 
L2  VpA/ 

61.70/El  \  1/2 
L2  VpA/ 


Figure  Cl.  First  three  modes 
and  frequencies  of  a  uniform 
cantilever  beam  based  on  ele¬ 
mentary  beam  theory  (after 
Warburton  1964) 


where 

0) 

o 

L 

E 

I 

P 

A 

Mode  1 


natural  frequency  of  beam  excluding  secondary  effects 
length  of  cantilever  beam  (185  ft) 
modulus  of  elasticity  (576,000  ksf) 

4 

moment  of  inertia  (42,667  ft  ) 

3 

weight  density  of  material  (0.150  k/ft  ) 

2 

effective  total  cross-sectional  area  (80  ft  ) 


w 


3.516 

2 

185 


'576, 000  x  42,667  \ 
k  0.15  x  80  / 


1/2 


4.649  cyc/sec 


*  References  cited  in  this  appendix  are  included  in  the  references  at  the  end 
of  the  main  text. 
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Mode  2 


0) 

o 


22.03  ( 576,000  x  42,667' 


1/2 


185 


2  V  0.15 


80 


29.13 


Mode  3 


U) 

o 


61.70  / 576,000  x  42,667\  1/2 
1852  V  0.15  x  80  ; 


81.59 


Timoshenko's  Theory 

2.  This  discussion  concerns  the  natural  frequency  of  a  cantilever  beam 
including  the  effect  of  transverse  shear  and  rotary  inertia  (Kruszewski  1949). 
The  effect  of  shear  lag  and  shear  deformation  of  the  web  is  to  increase  the 
flexibility  of  the  beam  because  of  the  additional  deflection  that  is  intro¬ 
duced.  The  effect  of  rotary  inertia  is  to  increase  the  dynamic  loading  on  the 
beam  because  of  the  additional  inertia  loading  due  to  the  rotational  acceler¬ 
ation  of  the  differential  elements  of  the  beam.  Considerable  lowering  of  the 
frequency  due  to  the  secondary  effects  is  obtained  for  the  higher-mode  num¬ 
bers,  as  shown  in  Figure  C2  (Kruszewski  1949) 


where 


u>o  =  natural  frequency  of  beam 

k  =  frequency  coefficient  where  shear  and  rotary  inertia  are  neglected 

D 

O 


m  =  mass  of  beam  per  unit  length  , 
k  =  coefficient  of  shear  rigidity  I  -jr 


A  =  shear  area 
s 


u  =  Poisson  ratio  (0.20) 

E 


G  =  shear  modul 


(f  AI  ■  66-7) 

io  (0.20) 

"s  [lu 


El 

A  G 
s 


=  0.21 


+  u) 


=  240,000  ksf 
1 


k^  =  coefficient  of  rotary  inertia  !jr 


=  0.12 


A  =  effective  total  cross-sectional  area  (80  ft  ) 
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(a)  FIRST  MODE 


(b)  SECOND  MODE 


(c)  THIRD  MODE 


Figure  C2.  Illustration  of  how  to  determine  the  ratio  and  natural 
frequencies  of  a  cantilever  beam  with  and  without  considering  shear 
and  rotary  inertia  (From  Kruszewski  1949) 


Mode  1 


—  =  0.885/.  a)  =  0.885(4.649)  =  4.11  cyc/sec 


Mode  2 


--  =  0.58  .'.u  =  0.58(29.13)  =  16.90 


Mode  3 


_  =  0.44  .-.w  =  0.44(81.59)  =  35.90 
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APPENDIX  D:  GTSTRUDL  INPUT  AND  OUTPUT  FILES  FOR 
GRAVITY  DAM  EXAMPLE  PROBLEM 


3TRUDL  '  RBRDAT 1 '  'RICHARD  B.  RUSSELL  DAM  NON-OVERFLOW  MONOLITH' 

♦ 

UNITS  FEET  KIPS 
$ 

*  FIRST  GENERATE  ALL  JOINTS  AS  HAVING  ZERO  COORDINATES  AND  THEN 
$  GENERATE  THE  CORNER  NODE  COORDINATES  FOR  ALL  ELEMENTS  IN  THE 

$  CHANGES  MODE.  GTSTRUDL  WILL  ASSUME  THAT  THE  MIDSIDE  NODES 
$  HAVE  COORDINATES  OF  ZERO. 

* 

GENE  135  JOI  ID  1  1  X  0  0 
CHANGES 

JOINT  COORDINATES 
1 

9  143.25 
29  3.84  46. 

37  104.92  46. 

85  11.92  143. 

93  40.25  143. 

99  11.92  160. 

107  28.92  160. 

113  11.92  170. 

121  28.92  170. 

127  11.92  185. 

135  28.92  185. 

$ 

GEN  B  1  9  37  29 

XD  4  P  EQ 

YD  2  P  EQ 

GEN  B  29  37  93  85 

XD  4  P  EQ 

YD  4  P  EQ 

GEN  B  85  93  107  99 
XD  4  P  EQ 
YD  1  P  EQ 

GEN  B  99  107  121  113 
XD  4  P  EQ 
YD  1  P  EQ 

GEN  B  113  121  135  127 
XD  4  P  EQ 
YD  1  P  EQ 
$ 

ADDITIONS 

$ 

TYPF  PT  AUF  CTppcc; 

GENERATE  4  ELEMENTS  ID  1  1  F  1  2  T  3  2  T  1 7  2  T  1 5  2  T  2  2  T  1 1  1  T  16  2  T  10  1 

REPEAT  8  ID  4  F  14 

STAT  SUPPORT  1  TO  9 

CONSTANTS 

E  576000.  ALL 

POISSON  0.20  ALL 

HEN  0.150  ALL 

EI.EM  PROP 

1  TO  36  TYPE  '  I PPQ '  THICK  1.0 
DAMPING  0.05  4 
l 

;  IN  ITS  KIPS  FEET  SECONDS  CYCLES 


SEED 


DUMP 


STORE  RESPONSE  SPECTRA  ACCELERATION  LIN  VS  PERIOD  LIN 
$  ACCELERATION  (FT/SEC**2)  VS  PERIOD  (SEC) 


DAMPING 

0 . 05 

FACTOR 

0.25 

31.78 

.0050 

33.42 

.0260 

51 . 23 

.0745 

59.70 

.  0805 

78.99 

.  1260 

81.69 

.  1635 

72.61 

.  2680 

68.68 

.  2899 

57.06 

.3755 

53.19 

.  4005 

40.22 

.5055 

37 . 22 

.5415 

27.05 

.  6980 

23.92 

.  7955 

17.87  1 

.  0905 

16.00 

1  .  1945 

9.98  1 

.6610 

8.82 

1 .8030 

END  OF  RESPONSE  SPECTRUM 
$ 

UNITS  KIPS  FEET  SECONDS 
INERTIA  OF  JOI  LUMPED 
INERTIA  OF  JOINTS  MASS 
$  HYDRODYNAMIC  'ADDED  MASS' 
1  TRANS  X  1.55 


10 

TRANS 

X 

3.02 

15 

TRANS 

X 

2.91 

24 

TRANS 

X 

2.79 

29 

TRANS 

X 

2.74 

38 

TRANS 

X 

2.68 

43 

TRANS 

X 

2.53 

52 

TRANS 

X 

2.37 

57 

TRANS 

X 

2.20 

66 

TRANS 

X 

2.01 

/I 

TRANS 

X 

1.81 

30 

TRANS 

X 

1 .58 

85 

TRAnS 

X 

1  .  14 

94 

TRANS 

X 

.76 

99 

TRANS 

X 

.46 

108  TRANS  X  .23 
113  TRANS  X  .06 

$ 

LOAD  1  'HYDROSTATIC  PRESSURES' 

ELEMENT  LOADS 

1  EDGE  FOR  EDO  4  GLO  VAR  VX  9.17  9.39  10.61 

5  EDGE  FOR  EDG  4  GLO  VAR  VX  7.74  8.4G  9.17 

9  EDGE  FOR  EDG  4  GLO  VAR  VX  6.22  6.98  7.74 

13  EDGE  FOR  EDG  4  GLO  VAR  VX  4.71  5.47  6.22 

17  EDGE  FOR  EDG  4  GLO  VAR  VX  3.20  3.96  4.71 

21  EDGE  FOR  EDG  4  GLO  VAR  VX  1.69  2.44  3.2 

25  EDGE  FOR  EDG  4  GLO  VAR  VX  .62  1.15  1.69 

29  EDGE  FOR  EDG  4  GLO  VAR  VX  0.0  .31  .62 

LOAD  2  'DEAD  LOAD' 

ELEMENT  LOADS 

1  TO  36  PODY  FORCES  GLOBAL  BY  0.150 
UNITS  LBS  INCHES 
STIFFNESS  ANALYSIS 

DYN  LOAD  3  'SEED  RESPONSE  SPECTRUM' 

SUPPORT  ACC 

TRANSLATION  X  FILE  'SEED' 

END  OF  DYN  LOAD 


2ND  OF  DYN  LOAD 

s 

E  IGENPROBLEM  PARAMETERS 
SOLVE  USING  SUBSPACE  ITERATION 
NUMBER  OF  MODES  4 

PERFORM  NO  STURM  SEQUENCE  CHECK  ORTHOGONALITY  CHECK 
TOLERANCE  EIGENVAL  l.E-4 
END 
$ 

DYN  ANAL  MODAL 

COMPUTE  DYNAMIC  STRESSES  MODAL  COMB  ALL 

PRINT  DYN  DATA 

OUTPUT  DECIMAL  4 

OUTPUT  BY  LOADING 

OUTPUT  FIELD  E 

CREATE  PSEUDO  STATIC  LOADING  4  'CQC  OF  LOADING  3'  AS  CQC  OF 
DELETIONS  ;  LOAD  3  ;  ADDITIONS 

LOADING  COMBINATION  5  'STATIC  LOAD  =  DEAD  LOAD  +  HYDROSTATIC 
COMBINE  5  1  1.0  2  1.0 

LOADING  COMBINATION  6  'STATIC  +  DYNAMIC' 

COMB  INE  6  5  1.0  4  1.0 
LIST  STRESSES 

CALCULATE  AVERAGE  STRESSES 
SAVE  DIRECT  ' R  B  R  S  A  V 1 ' 

F  INISH 
END  OF  FILE 


AAAAAAAAAAAAAAAAAAAAAAAAA 
A  EIGEN-SOLUTION  CHECKS  A 
AAAAAAAAAAAAAAAAAAAAAAAAA 
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MODE-- 

- EIGENVALUE - 

( (RAD/SEC) AA2) 

- FREQUENCY - 

( RAD/SEC  > 

- FREQUENCY - 

(CYC/SEC) 

] 

1.441 404D+03 

3 . 7  9  6  5  8  3  D  +  0  1 

6 . 042449D+00 

1  . 

1 

7 . 008182 D +03 

3 . 371 488D+0 1 

l . 3 32364D+01 

7  . 

1 . 745866 D +04 

1 .321312 D +02 

2 . 1 02934D+0 1 

4  . 

•1 
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2. 237143 D +01 

4  . 
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LOADING  3 


-PERIOD - 

l SEC/CYC) 

654958D-0 1 
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APPENDIX  E:  COLE/CHEEK  COMPUTER  PROGRAM  INPUT  AND 
OUTPUT  FILES  FOR  GRAVITY  DAM  EXAMPLE  PROBLEM 


«OLD  G2DGN I/SDFDAM1 ,R 
tFRN 


DO  VOU  UISH  TO  SEE  I NFORM AT  ION  FILE  ? 

VES  <Y>  OF  NO  <CR) 

*N 

NCOL,  B0  (-NROU  IF  SM  -0.),  SM  7 
■-4,9,0 

INPUT  DAM  GEOMETERV  IN  FT  AND 

ELEUATIONS  RELATIUE  TO  ANY  DATUH 
ELEU.  OF  BASE  7 

•  0.0 

SLOPE  s  RUN  TO  RISE  RATIO 
UPSTREAM  SLOPE  ? 

-0.9833 

BREAK  ELEU.  OF  UPSTREAM  SLOPE  ? 

•143. 

DAM  CREST  ELEU.  ? 

■18S  . 

CREST  WIDTH  ? 

•  17  . 

DOUNSTREAM  SLOPE  7 
•0 .67253 

BREAK  ELEU.  OF  DOUNSTREAM  SLOPE  ? 

-1  60. 

RADIUS  OF  TRANSITION  7 

•0 . 

SUFPLY  NROU+ 1  UALUES  FOR  YH 

•  0 •  S3.  46.  70.25  94. S  118.75  143.  160.  170. 


E3 


185  . 


INITIAL  PHASE  -  PROGNESS  INDICATOR 

. >1/5  * 

. >2/5  * 

. >3/5  * 

.  .  .  . . >4/5  * 


MODULUS  OF  ELASTIC  ITV  *  E  (MILLION  PSI)  OR  F'C  (PSI) 
■  4  . 

UNIT  UEIGHT  OF  CONCRETE  t  GANC  ( LBF/CF  ) 

-  150  . 

RESERVOIR  ELEVATION  i  HUATER  (FT) 

•170. 


NATURAL  PERIOD  OF  DAN  •  0.12950  SECONDS 

NATURAL  PERIOD  OF  DAN  ♦  UATER  •  0.15983  SECONDS  (  6.28  HTZ ) 


DO  VOU  DISH  TO  USE  SEED'S  (MEAN) 

DESIGN  SPECTRUM  (SX)  FOR  ROCK  SITES  7  (CR) 
SUPPLV  VOUR  OUN  SPECTRUM  FILE  7  (1) 
SUPPLY  SPECTRAL  ACCELERATION  VALUE  7  (2) 

(ZERO  VALUE  FOR  STATIC  STRESSES  ONLY) 

■ 

PEAK  GROUND  ACCELERATION  ? 

-0.25 


SOLUTION  PHASE  -  PROGRESS  INDICATOR 

. >1/5  * 

. >2/5  f 

. .  t 


E4 


NROU  -  9  NCOL  •  4  NEQ  -  90  MB AND  •  14  NB  LOCK  •  i 

. >4/5  * 


1509  HRS.,  14  MAY  1987 


DAM  HEIGHT  (FT)  1 
POOL  HEIGHT  (FT)  t 
MODULUS  (PSI  *  10**6)  l 
SPECTRAL  ACCELERATION  (G'S)  t 


185.0 

170.0 

4.00 

0.61 


PRINCIPAL  STRESSES  AS  A  FRACTION  F'C«  4353. 
WHERE  E  •  33.  *  ( GAflC  **  1.5)  *  SQRT(F'C) 


STRESSES  INCLUDE  GRAUITY  AND  HYDROSTATIC  LOADS 


UPSTREAM 

PRINCIPAL 

DAM  HEIGHT 

DOWNSTREAM 

PRINCIPA 

STRESSES  (XF'C) 

(FT) 

STRESSES 

(XF'C) 

36.  ( 

IX  ) 

177.5 

19. 

( 

0X  ) 

122.  ( 

3X  ) 

16S.0 

165. 

( 

4  X  ) 

268 .  ( 

6  X  ) 

151.5 

295. 

( 

7X  ) 

328.  ( 

8  X  ) 

130.9 

270. 

( 

6X  ) 

329.  ( 

8  X  ) 

106.6 

255. 

( 

6X  ) 

333.  ( 

8X  ) 

82.4 

237. 

( 

5X  ) 

344.  < 

8X  ) 

58.1 

208. 

( 

5X  ) 

379.  ( 

9X  ) 

34.5 

163. 

( 

4X  ) 

494.  ( 

11X) 

11.5 

114. 
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